Three-dimensional inverse energy cascade induced by vortex reconnections 
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A recent study of homogeneous isotropic turbulence by Biferale, Musacchio and Toschi has deter- 
mined that a three-dimensional inverse energy cascade is possible if the nonlinearity of the Navier- 
Stokes equation is restricted in Fourier space to helical modes of the same sign. In low-temperature 
superfluid helium, viscosity is zero, vorticity takes the form of discrete, thin vortex filaments of fixed 
circulation, and turbulence is a tangle of such filaments. We exploit the simpler nature of quantum 
vorticity to show that the three-dimensional inverse energy cascade can arise from reconnections of 
vortex loops of the same polarity which shift energy from small length scales to large length scales, 
in analogy to what was envisaged by Biferale and collaborators in classical Navier-Stokes turbu- 
lence. We discuss superfluid turbulence experiments and the observed generation of the classical 
Kolmogorov energy spectrum in view of this finding. 



PACS numbers: 67.25.dk, 47.37.+q 47.27.-i 



2 



The phenomenology of three-dimensional turbulence is based on Richardson's idea [l?} of the (forward) turbulent 
cascade: kinetic energy, injected externally, feeds large eddies, which interact, become stretched and deformed into 
smaller and smaller eddies, until, at sufficiently small length scales, viscous forces dissipate this energy into heat. A 
reversed flux of energy, from the small scales to the large scales, is observed in two-dimensional turbulence fl2l I22I]. 
Such inverse cascade is rare in three-dimensional turbulence, but can be observed in the presence of strong anisotropy 
M, E! HI , or when rotation [H, E 

[2l| or confinement in layers[24j make the flow almost two-dimensional. 



Recently, Biferale et al. [6[ numerically induced the three-dimensional inverse energy cascade by artificially restrict- 
ing the nonlincarity of the governing Navicr-Stokcs equation to the interaction of Fourier modes of the same helical 
sign. Their result clarifies previous observations and shows that, in principle, all three-dimensional turbulent flows 
contain nonlinearities which may lead to an inverse cascade: to make the effect apparent one has to break the mirror 
symmetry of the interactions. 

The aim of this report is to exploit the unusual properties of supcrfluid helium and consider the three-dimensional 
inverse energy cascade in a quantum fluid setting in view of the finding of Biferale et al. @ . At temperatures below 
1 K, liquid helium ( 4 He) is a pure supcrfluid. Unlike ordinary fluids (in which vorticity is a continuous field), the 
superfluid's rotational motion is constrained by quantum mechanics to discrete vortex lines of fixed circulation and 



atomic thickness. Vortex lines reconnect if they come sufficiently close to each other [111, [14J. In a recent experiment, 
Walmsley & Golov used high voltage tips to inject vortex rings into 4 Hc at very low temperatures, creating 
turbulence and monitoring its decay. Depending on the duration of the injection time, they found two decay regimes, 
characterized by distinct power laws. In a related experiment, Bradley et al. Q used a vibrating grid to inject vortex 
loops into the isotope 3 Hc-B (a fcrmionic supcrfluid), and observed the same power laws seen by Walmsley & Golov. 
We shall demonstrate that, in conditions typical of supcrfluid turbulence experiments, reconncctions of vortex loops 
with the same polarity induce a three-dimensional inverse energy cascade, thus providing Bifcrale's spectral argument 
with a simple, striking example in real physical space. 

We model vortex lines as space curves s(£,i), where £ is arc length and t is time, which move according to the 
Biot-Savart law [H|: 

ds n f (s — r) 

7J7 = -T- f T FT xdr ' 1 

dt An J c |s - r| 3 w 

where k = 9.97 x 10~ 4 cm 2 /s is the quantum of circulation in 4 He, and the line integral extends over the entire vortex 
configuration C. Our calculations are performed in a periodic cube of size D = 1 cm. The numerical techniques to 
discretize vortex lines into a variable number of points N held at minimum separation 8/2 = 0.01 cm, compute the 
time evolution, de-singularize the Biot-Savart integrals, evaluating them using a tree- method (with critical opening 
angle 0.4), and the algorithm for vortex reconnections, are all described in the literature fllfl fl9|. 

The initial condition at t = consists of vortex rings of random size taken from a distribution with mean radius 
R = 0.1 cm and standard deviation a = 0.047 cm; the rings, randomly placed in the computational box, are 
initially parallel to the xy-plane and travel in the positive z-direction. No further energy is injected into the vortex 
configuration. During the time evolution, the total energy slowly decreases because of numerical dissipation: the 
finite discretization along the vortex lines, controlled by the parameter 8, damps high frequency Kelvin waves, and 
the vortex reconnection algorithm, again controlled by S, removes vortex length. These numerical effects play a 
dissipative role similar to phonon emission in the real physical system (kinetic energy is turned into sound energy 
which is radiated away). In order to highlight the three-dimensional inverse energy cascade, the initial condition 
used here is simpler than the initial condition that we used to model [f| the experiment of Walmsley & Golov 23 1 



as realistically as possible (vortex rings were injected continually for a given period of time in the form of a narrow 
beam). 

We present the numerical results. Fig. [1] shows that, starting from the initial condition, the vortex rings interact, 
reconnect and quickly form a vortex tangle. The time evolution of the vortex line density L is shown in Fig. [5] We 
analyze the build-up of the energy spectrum E(k), which is defined by 



_ J -v 2 dV = J E(k)dk, (2) 

where v is the velocity field generated by the vortex lines, k = |k| is the magnitude of the three-dimensional wavenum- 
ber and V = D 3 is volume. Fig. [3] shows that at t = the energy is concentrated at k w 2ir/2R ss 30 cm -1 
corresponding to the initial vortex rings. Our main finding is that, during the evolution, the energy spectrum E(k) 
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FIG. 1. Snapshots of the vortex tangle at t — (left), t = 10 s (middle) and t — 20 s (right) 
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FIG. 2. Vortex line density L (cm 2 ) vs time f (s). 

increases at small wavcnumbers k < 30 cm" 1 and decreases at large wavenumbers k > 30 cm -1 . This transfer of 
energy from small scales to large scales (inverse energy cascade) is also visible in Fig. 3] where we plot the energy flux 
e(k) = dE(k')/dt dk' where k D = 2n/D. 

We argue that this inverse energy cascade arises from reconncctions of vortex lines. The interaction of two vortex 
lines occurs either indirectly via the velocity fields which they generate (the Biot-Savart law), or directly via vortex 
reconnections. To isolate the role played by vortex reconnections we repeat the calculation, replacing the Biot-Savart 
law with its Local Induction Approximation (LIA) [1, [l6| : 



^-ln(Ji/a„) S 'x S ", (3) 
at 4-7T 

where a prime denotes derivative with respect to arc length, R = l/|s"| is the local radius of curvature and do ~ 10 _8 cm 
is the vortex core radius in 4 He. Under LIA, the vortex lines ignore each others' velocity fields and interact only when 
they reconnect. Fig. [5] shows that the inverse energy cascade is still present under LIA. 

The following simple geometrical consideration captures the role of vortex reconnections on the scale-to-scale energy 
transfer. Energy and speed of a vortex loop of size R are roughly proportional to R and 1/R respectively. A head-one 
collisions of two vortex loops forms two loops of approximately the same size (see Fig. |B] top). A collision from behind 
of two vortex loops travelling in the same direction (see Fig. [6] bottom) forms a larger vortex loop which contains 
most of the energy (and which, being slow, is likely to be hit by other loops, thus becoming entangled), and a smaller, 
less energetic loop (which moves away and is likely to be absorbed by the walls). Clearly, the first type of collisions 
tends to leave length scales unaltered, and the second type tends to shift energy to larger length scales. 

When modelling the experiment of Walmsley & Golov [lH @, we found that a short injection time was followed 
by a decay of the form L ~ t _1 , and a long injection time was followed by a decay of the form L ~ t~ 3 / 2 , as observed 
in the experiment. It has been argued that the first decay is typical of random superfluid turbulence, such that the 
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FIG. 3. (Color online) Energy spectrum E(k) (arbitrary units) vs wavenumber k (cm' 1 ) at t — (dashed blue line), t = 10 s 
(dot-dashed red line) and t — 20 s (solid black line). Note the progressive increase of kinetic energy at small wavenumbers 
k < 30 cm -1 and the decrease at large wavenumbers. The energy spectra are computed from a 512 2 Cartesian mesh in the 
xz-plane. 




FIG. 4. Time averaged energy flux, (e) (arbitrary units) vs wavenumber (cm 1 ). 

only length scale present is the average intervortex distance I ~ L" 1 / 2 , and the second decay of more structured, 
quasi-classical turbulence [2(J ■ 

We now examine what happens if vortex loops are injected continually (by a vibrating grid for example), rather 
than given as initial condition and allowed to decay. Starting from an empty computational domain, we inject vortex 
rings with size drawn from the same normal distribution used earlier (R = 0.1 cm, a = 0.047 cm) with a frequency of 
125Hz. The rings are injected in the yz-plane and travel in the positive x direction; their evolution is computed with 
the full Biot-Savart law. Figure [7] shows snapshots of the vortex tangle as the simulation progresses. We track the 
evolution of the energy spectrum, and find that the inverse cascade is very capable of transferring sufficient energy 
to low k that a spectrum consistent with the Kolmogorov fc -5 / 3 scaling is built up, see Fig. [5J We remark that, 
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FIG. 5. (Color online) Energy spectrum E(k) (arbitrary units) vs wavenumber k (cm^ 1 ) at t — (dashed blue line), t = 10 s 
(dot-dashed red line) and t = 20 s (solid black line) for the calculation performed with LIA. Note again the progressive increase 
of energy at small wavenumbers. 




FIG. 6. Schematic head-on collision of vortex rings travelling in opposite direction (red, top) and collision from behind of 
vortex rings travelling in the same direction (blue, bottom), and resulting vortex configuration after the vortex reconnection. 



in contrast to our early work [5| which modelled the conditions of the experiment of Walmsley & Golov [23], the 
injection is not in the form of a narrow beam from a single point (which may induce a large scale motion by itself), 
but is uniform in the j/z-plane: the build-up of energy at small k is thus due only to the inverse cascade. Once the 
Kolmogorov spectrum is formed, if the injection of rings is stopped, the vortex line density decays as L ~ t~ 3 / 2 , as 
observed in the experiments 0, H^j and in numerical calculations [H[ . 

In conclusion, superfiuid turbulence experiments such as those of Walmsley & Golov [23| and of Bradley et al. 
involve injecting a beam of vortex loops into the experimental cell, creating an initial anisotropy. This set-up 
creates a condition which favours the interaction and reconnection of vortex loops of the same polarity (which tends 
to transfer energy to larger length scales) and suppresses the head-on interaction of loops of the opposite polarity 
(which tends to leave the length scales unchanged). In physical space, this selection mechanism is analogue to the 
decimation procedure applied by Biferale et al. [6[ in Fourier space to the Navier-Stokes equation, which suppressed 
the interaction of helical modes with different sign, and induces the inverse energy cascade. 

We acknowledge fruitful discussions with W.F. Vinen, J. Laurie and S. Nazarenko, and the financial support of the 
Lever hulme Trust and the EPSRC. 
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FIG. 8. (Color online) Energy spectrum E(k) (arbitrary units) vs wavenumber k (cm x ) at t = 0.75 s, L = 43.9 cm 2 , (dashed 
blue line) and t = 15 s, L = 159.6 cm -2 , (solid black line) corresponding to a continual injection of vortex rings, see Fig. [7] 
Note the formation of the Kolmogorov fc~ 5 ^ 3 spectrum, indicated by the dot-dashed red line. 
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